Steady state entanglement of two superconducting qubits engineered by dissipation 
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We present a scheme for dissipative preparation of an entangled steady state of two superconducting qubits 
in a circuit QED setup. Combining resonator photon loss, a dissipative process already present in the setup, 
with an effective two-photon microwave drive, we engineer an effective decay mechanism which prepares a 
maximally entangled state of the two qubits. This state is then maintained as the steady state of the driven, 
dissipative evolution. The performance of the dissipative state preparation protocol is studied analytically and 
verified numerically. In view of the experimental implementation of the presented scheme we investigate the 
effects of potential experimental imperfections and show that our scheme is robust to small deviations in the 
parameters. The promising results of our study thus open a route for the demonstration of an entangled steady 
state in circuit QED. 
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I. INTRODUCTION 

One of the most peculiar properties a physical system can 
exhibit is quantum-mechanical entanglement 1 1]. From a fun- 
damental perspective, entanglement is a non-classical effect 
which is indispensable for the understanding of fundamental 
quantum physics. From a technological perspective, entan- 
glement is useful for enhanced measurement techniques and 
is an important element in quantum information processing 
and quantum communication \2\. For the past two decades 
great effort has been invested into the generation and inves- 
tigation of entangled states. Inspired by the circuit model of 
quantum computation entanglement has predominantly been 
investigated by means of coherent interactions, i.e. by apply- 
ing sequences of unitary gates. Today, there is a large num- 
ber of physical systems where entanglement has been demon- 
strated and which are considered suitable for the realization 
of advanced quantum information protocols. Out of these, su- 
perconducting systems |3| have proven to be good candidates 
for the realization of quantum algorithms involving many gate 
operations [^-§\. Despite impressive reductions of the deco- 
herence Q in superconducting systems II8 UT2JI . any state other 
than the ground state will deteriorate over time. As a conse- 
quence, any entangled state produced by gate operations will, 
at most, be meta-stable. 

Over the past few years, however, an alternative approach 
of dissipative state engineering and dissipative quantum com- 
puting lfT3l [141 has emerged and gained increasing attention. 
As opposed to unitary quantum computing, where decoher- 
ence and dissipation act detrimentally on the state preparation 
process and on the prepared state, the central idea here is to 
prepare non-trivial quantum states relevant for quantum in- 
formation, simulation ifTSl lT&l. memories ifTTl . or communi- 
cation [18 1 by the dissipative interaction of the system with 
its envkonment. As a result, a stationary state is obtained 
that is resilient to the sources of decoherence by which it has 
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been produced. This idea has been studied theoretically for a 
number of quantum optical and solid state systems, in paiticu- 
lar cavity QED 1 19-25], atomic ensembles f26'-'28l, ion traps 
1 16, 29-31], plasmonic systems |32-34| and optical lattices 
IIT5II35I . The first experimental demonstrations were achieved 
in atomic ensembles |36| and ion traps |31 1. 

Recently, proposals utilizing dissipation have also been de- 
veloped for superconducting systems Il37l - l40l . These are ei- 
ther aiming for the preparation of steady states which exhibit 
a certain degree of entanglement |37|, squeezed states |[38l , or 
entangled states in the presence of ultrastrong coupling |39|. 
Most recently, cavity cooling into arbitrary single qubit states 
has been studied both theoretically and experimentally |40|. 

In this work, we present a scheme for the dissipative prepa- 
ration of a maximally entangled steady state of two super- 
conducting systems which are strongly coupled to a mode 
of a common resonator. Our scheme is specifically designed 
to exploit the level structure of typical transmon qubits I.41J , 
which constitute a weakly anharmonic oscillator. The scheme 
is, however, not particularly restricted to transmons, but one 
could also consider phase qubits |42| coupled to a resonator. 
As opposed to previous dissipative entangling schemes for 
cavity QED 12311251 . where the cavity couples only to a sin- 
gle transition in a A configuration, one cannot restrict the cav- 
ity coupling to a single transition in the transmons due to the 
weak anharmonicity. In this article we show that not only is 
it possible to make a scheme suitable for transmon qubits, the 
coupling to several transitions in the transmon is in fact an 
advantage since it automatically creates a certain reshuffling 
of undesired states, which would otherwise have to be imple- 
mented by additional driving fields. As a consequence, our 
scheme works ideally for realistic transmon qubits where the 
anharmonicity is rather low, which is a requirement for having 
good coherence properties pT) . 

In the setup, detailed in Sec. |llj we exploit the strong dissi- 
pation through the resonator. In addition, we utilize a coherent 
two-photon drive of a dipole-forbidden transition with a two- 
tone microwave field, similar to Refs. |10 43 1. Using these 
ingredients, we engineer an effective resonator loss process 
which deterministically prepares a maximally entangled sin- 



glet state, as described in Sec. Ill This is achieved by driving 
the system into a particular excited state which rapidly decays 
to the singlet. This process is resonant by construction, as 
opposed to any processes that take the system out of the sin- 
glet which are off -resonant and therefore strongly suppressed. 
Given that the target state is produced by a time-independent 
loss process and continuous wave fields, it is a steady state of 
the dissipative evolution. 

In Sec. IV we investigate the performance of our scheme, 
both analytically, to derive benchmarks for the protocol, and 
numerically, to verify the mechanisms that underlie the pre- 
sented dissipative state preparation scheme. In view of an ex- 
perimental realization of our proposal, we also study effects 
of imperfections of the system. Our results show that a max- 
imally entangled state of two superconducting qubits can be 
prepared rapidly and with a high fidelity for realistic imperfec- 
tions in current circuit QED experiments. The proposal devel- 
oped here thus opens a route for the dissipative preparation of 
maximally entangled states of superconducting systems using 
existing technology. 



II. SETUP: COHERENT AND DISSIPATIVE 
INTERACTIONS OF TWO COUPLED TRANSMONS 

For our study we consider two superconducting transmons 
II4TI coupled to a common resonator in a circuit QED setup. 
The coherent dynamics of the system is described by a Hamil- 
tonian H = iJfrcc + ^cav + ^d- The energy levels are illus- 
trated in Fig. [Tja) and described by the free Hamiltonian 
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J2 (2c.-2A)|2),(2|+^|l),(l|, (1) 



with levels | k) of transmon j and the resonator mode a. Here, 
u! denotes the level spacing of the two lower levels and A the 
anharmonicity, with h = 1. In our analytical discussion we 
will focus on the first three levels of the transmons; further 
levels are included in our numerical simulations. 

The transitions of the transmons, |0) O |1) and |1) o |2), 
are coupled by the coherent interactions shown in Fig. [lib). 
They are described by a Hamiltonian -ffcav + ^d- Here, TJcav 
represents the coupling of the resonator to the transitions of 
the transmons, 

ffcav = Y. 5«^ (|0)j(l| + ^|l)j(2|) +H.C., (2) 

3 = 1,2 

with a coupling constant g, and a factor of \/2 for the matrix 
element of the upper transition. The coherent drive 



J = l:2 



x(|1),(0| + n/2|2),(1|) +iJ.c. 



(3) 



contains several microwave fields which couple the transi- 
tions |0) <-^ |1) and |1) i^ |2). We assume that the drive 
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FIG. 1: (Color) Setup. The internal levels of two transmons (a) are 
coupled by coherent interactions (b) to mimic the A system in (c). 
Two microwave fields fii/2 provide virtual couplings of the transi- 
tions |0) -O- |1) and |1) <->■ |2) (b) which combine to an effective 
two-photon drive Jloff of the transition |0) -f-)- |2). The transmon- 
resonator coupling (g) is resonant with the upper transition and de- 
tuned by (5i — 5c from the lower transition. Spontaneous emission (7) 
and resonator photon loss (k) are present as decoherence processes. 
The detunings are defined in the text. 



with f22 exhibits an opposite phase for each of the two trans- 
mons, whereas the phase of fJi is identical for both transmons. 
As we will see, this choice of phases allows us to break the 
symmetry of the system and thereby drive certain transitions 
which play an important role in our proposal. 

We choose the frequencies of the two fields in such a way 
that they combine to an effective two-photon drive of the tran- 
sition |0) O 1 2) with a coupling constant of fioff that will be 
derived in Sec. |IIIA In doing so, we render the couplings 
of the system resembling the A system shown in Fig. [T] c), 
with (meta-) stable lower levels |0) and |1) and an "excited" 



level 1 2) for each of the transmons. "Excitation" from |0) to 
1 2) is then accomplished by the two-photon drive with fioff. 
For most of this paper, we will assume that the resonator cou- 
pling is resonant with the transition |1) f-> |2), while being 
somewhat detuned from the lower transition |0) -^ |1). 

In the following, we will avoid the fast dynamics in the 
drive by changing into a frame rotating with a Hamiltonian 
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FIG. 2: (Color), (a)-(b) Dissipative state preparation mechanisms, (c) loss mechanisms and (d) effective lower-level decay processes, (a) 
Effective resonator decay from |00) into |S) involves coherent coupling to |So). jSo) and |S)|1) are strongly coupled (\/25) so that these states 
hybridize and form dressed states | S± ) (shown here for a choice of 5c = 82 — 5i, by which the resonator is resonant with the upper transition). 
By setting 82 — \/2<7 the driving from jOO) is resonant with the lower dressed state |S_). Population from |00) is thus rapidly excited and 
decays into |S) via the effective engineered resonator decay k+. (b) The population of the states |11) and jT) is shuffled to |00)|0) by the 
resonator coupling g and successive resonator decay at an effective rate of Kcff . (c) The two-photon drive also causes an undesired coupling 
of |S) to an excited state |Ti). |Ti), in turn, couples to a number of (resonator-) excited states which form dressed states at different energies 
(indicated) and eventually decay to other states. These can generally be made off-resonant with the drive from |S) by an appropriate choice 
of the resonator and microwave detunings so that the effective resonator decay k- from jS) is suppressed. In addition, |S) decays into |00) 
through qubit decoherence 7. The effective decay processes of the lower levels are summarized in (d). 



where w = | (wi + W2) is the mean frequency of the classical 
driving fields. Applying a unitary Z^ = exp[iiJiot^] we obtain 
a transformed Hamiltonian H = UHW + iLiW = H, 



H^ 



H^ in a frame rotating with Tfic 



free ' 

The transformed 



free Hamiltonian can be expressed as 



h',,,, = Wa+ V 6i\l),{l\+S2\2),{2\, (5) 



3=1,2 



where Si = uj — ui, 62 = 2(cj — w) — 2 A, and (5c = Wc — a) 
denote the energies of the transmons and the resonator in the 
rotating frame. Furthermore, we obtain the interaction Hamil- 
tonians H^^^ — if^av for 'h^ transmon-resonator coupling and 
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for the drive. With this choice of the reference frame rotating 
with the mean frequency, we find for the detunings of the 
microwave fields A1/2 = w — wi/2 = ±ij (w2 — wi). 

In addition to the coherent dynamics discussed so far, the 
system also exhibits dissipative couplings, which is essential 
for the dissipative state preparation mechanisms we would 
like to engineer. The dissipative dynamics of the open sys- 
tem is determined by its coupling to a bath and the properties 
of the bath. Assuming the bath to be Markovian, the system 
dynamics is governed by a master equation of Lindblad form 
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with one Lindblad operator Lk for each physical decay pro- 
cess present in the system. As illustrated in Fig. [T] a), we 
assume that transmon j undergoes spontaneous decay which 
in the transmon regime can be described by 



^7ij=V7|0),(l| 



L^2,j = V27|l>,(2| 



(8) 
(9) 



We numerically verify that, with reasonable numbers 17-91, 
the effect of additional, pure dephasing is negligible. The pho- 
ton loss out of the resonator is described by 



L« 



'KO, 



(10) 



where k is the photon loss rate. 

Due to our choice of the couplings similar to a A config- 
uration, most of the dynamics will happen in the two lower 
levels. To describe them we choose a two-atom basis with 



triplet states |00) = |0)i|0)2, |11), |T) 



and the singlet state |S) 



V2 



(|oi> + |io)). 
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(|01) - 1 10)) as the desired en- 



tangled steady state. For the detailed discussion of the engi- 
neered decay processes, we also introduce the excited atomic 

states |To) = ^ (|02) + |20)), |So) = ^ (|02) - |20)), 

|Ti) = ^ (|12) + |21)) and |Si) = ^ (|12) - |21)). The 
presence of resonator excitations is indicated by a second ket 
vector, e.g. |00) 1 1). We use this notation to explain the mech- 
anisms of our scheme in Sec. HHlbelow. 



III. MECHANISMS FOR DISSIPATIVE PREPARATION OF 
THE MAXIMALLY ENTANGLED SINGLET STATE 



(7) In this section we will show how to engineer effective decay 

processes which prepare a steady state close to the maximally 



entangled singlet state |S). For now, we will focus our discus- 
sion on the physical mechanisms behind the effective decay 



processes, while Sec. Ill A and IIIB will deal with the deriva 



tion of quantitative expressions for the effective operators and 
rates. 

The mechanism of our scheme is illustrated in Fig. l2] a). 
The working principle is based on a strong asymmetry be- 
tween the decay into the singlet state | S) at a rate of k^ and the 
loss process from the singlet state at a rate of k_. In the fol- 
lowing we will discuss the physical mechanism for the prepa- 
ration of |S). 

In the previous section we have introduced a coherent driv- 
ing if (J. The purpose of it is to drive a two-photon transition 
|0) O 1 2). For now, we will assume that we have a coherent 
drive of |0) 4-> |2) with a coupling constant of flcs and defer 
the derivation to later Due to the opposite phase of fl2 on the 
two transmons, this drive then couples |00) to an excited state 
|So) with a detuning of 82, as can be seen from Fig. |2]a). |So) 
is in turn coupled to |S)|1) by the resonator coupling i/cav 
From here, |S) |1) decays into |S) via resonator decay at a rate 
of K. These processes combine to an effective resonator decay 
process from |00) into |S) with a rate of k+. 

In order to engineer this process to be as strong as possible 
we have to fulfill two requirements: First, we need to make 
sure that the coupling of the transmon-excited state |So) to the 
resonator-excited state |S)|1) is close to resonance, given that 
only the latter can decay to |S) through resonator photon loss. 
To this end we set the resonator into or close to resonance 
with the upper transition of the transmons, |2) |1)|1). This 
is reached by choosing ojc = w — 2A (S^ ^ 62 — 61), and 
results in an equal energy of |So) and |S) |1), as shown in Fig. 
l2]a). The two states hybridize and form dressed states 



|S±) 



1 



V2 



(|So)±|S)|l)), 



(11) 



located at frequencies of 2a; — 2 A ± \/2<7 (or $2 ± V^g)- 

The second requirement is that the two-photon drive from 
from |00) is resonant with one of the dressed states in Eq. (fTTli. 
Choosing a detuning of 62 = \/2g, we tune the drive into res- 
onance with the transition from |00) to |S_). Population from 
|00) is then rapidly excited to |S_), which, through its con- 
tribution from |S)|1), decays into |S). For a strong resonant 
drive, the resulting effective decay process is only limited by 
the line width | of |S_) as the state mediating it. Thus, the 
dissipative preparation mechanism of the singlet and its rate 
K+ can be engineered to be rather large. 

Loss from the singlet can occur through the couplings of |S) 
to any excited state other than |So) by the available microwave 
fields, e.g. to |Ti) by fioff- As indicated in Fig. |2]c), these 
excited states are coupled to a number of other, in particular 
resonator-excited states, for instance |Ti) couples to |11)|1), 
|To)|l), |T)|2) and |00)|3). Consequently, this estabHshes a 
loss channel from |S) through effective resonator decay, e.g. 
into 1 11), which causes losses at a rate k_ from the desired 
steady state |S). Fortunately, the photon-number dependent 
coupling strength between transmons and resonator provides 
us with a non-equidistant spectrum which consequently makes 



it possible to have the two-photon drive resonant with the tran- 
sition from 1 00) to |S_) while keeping it off -resonant with the 
transitions from |S) to other hybridized excited states. In this 
way, loss processes from the singlet are suppressed by their 
detunings. 

In order to reach |S) independently from the initial state 
and to maintain it as the steady state, an additional mecha- 
nism is required to transfer population from lower states other 
than |00), i.e. from |T) and |11), to |S). We have assumed 
that the resonator is resonant with the upper transition. Due to 
the anharmonicity, the resonator is off -resonant with the lower 
transition. For reasonable anharmonicities the off-resonant 
coupling is still sufficient to allow a reshuffling of population 
from 1 11) to |T) and |00). As is shown in Fig. |2]b), this reshuf- 
fling process involves the resonator coupling of the lower tran- 
sition (\/25), e.g. |T) <-> |00) 1 1), and decay of a resonator ex- 
citation at a rate of k. It can be seen as an effective decay pro- 
cess with a decay rate KcS — 2K.g^/[2,g^ + (Si — Sc)'^ + k^/4]. 
This expression contains both limiting cases, where one can 
either eliminate the resonator-excited states, or where the 
states can be seen as a dressed states with resonator-excited 
states, for instance the triplet states 



V2 



(|T)±|00)|1)), 



(12) 



which decay towards |00) at rates ex k. Ideally, the reshuffling 
mechanism rapidly transfers population of the triplet states 
and losses from the singlet (into |11)) to |00), from where 
they decay into |S) by the dissipative preparation mechanism 
discussed above. The fastest reshuffling is reached by tun- 
ing the resonator into resonance with the lower transition, i.e. 
(5c = Si. This choice is, however, different from the above 
choice of Sc = S2 — ^1 which optimizes the dissipative state 
preparation process. With this choice of the resonator fre- 
quency we get Kcff — '2'Kg'^/[2g^ + 4A^ + k^/4], from which 
we see that the reshuffling works best for small anharmonicity 
A. For larger A the process becomes less effective. Having 
both processes, state preparation and reshuffling, simultane- 
ously active might therefore seem problematic for large anhar- 
monicities. However, as we shall see below, the scheme can 
still be effective for large A if we allow for an longer time for 
the reshuffling. Furthermore, as we will also address below, 
the two requirements for Sc above are far less critical than the 
resonant set-up of the two-photon drive. Consequently, both 
processes, the dissipative state preparation and the reshuffling, 
can be effective at the same time over a wide parameter range. 



as we will numerically demonstrate in Sec. IV 



In addition to effective resonator decay, qubit decoherence 
present in the system can cause loss from the singlet indepen- 
dent of the drives. Most notably, it can cause a loss from |S) 
into |00), as shown in Fig. |2]c). We will see below that its 
effect is much less than the losses through effective resonator 
decay for realistic parameters. The presented mechanisms are 
summarized in Fig. l2]d): On the left hand side we see the 
reshuffling mechanisms enabled by the resonator coupling to 
the lower transition, represented by Kcb, and on the right hand 
side the state preparation (k+) and loss (k_) mechanisms af- 
fecting the singlet state, as well as the decay by from | S) by 



qubit decoherence at a rate of 7. 

To sum up this section, we have identified suitable mecha- 
nisms for the dissipative preparation of the singlet state and 
discussed the physical effects behind them. In the follow- 
ing two sections we will analytically derive the couplings and 
rates for the effective coherent and dissipative processes in 
our scheme. Based on these, we derive benchmarks for the 



performance of the scheme in Sec. IV 



A. Effective coherent driving of the dipole-forbldden 
transition 0) <^ 2) by a two-photon process 

The implementation of the dissipative state preparation 
scheme discussed above requires a coherent coupling of the 
transition |0) ^ |2). Since this transition is dipole-forbidden, 
such a coupling cannot be accomplished in a single step. One 
way to overcome this is to use a two-photon process, achieved 
by the combination of two individual fields. In H^ we have 
chosen two such fields, fli and ^2- As we will derive in the 
following, these provide complementary virtual single-photon 
excitations which form the desired coupling. 

In the following, we will apply the effective operator for- 
malism presented in Ref. |44| to obtain a simple effec- 
tive Hamiltonian for a single transmon + two-photon drive. 
Here, we separate the Hamiltonian into a perturbative part 
V{t) = H^, which contains the fields, and a perturbed part 



Ho — Hf 



Sco) a. (Note that the derivation below is for a 
single transmon only. With this in mind, the reuse of Hamil- 
tonian definitions should not cause any confusion.) While in 
Ref. Il44l only effective processes with an initial excitation are 
considered, here we also allow for an initial deexcitation. We 
therefore set up the effective Hamiltonian (cf. Ref. |i44J) as 
H,s ^ Hip + Hi-^Kiih 



/=lfe=0 



H 



(fc,/,±) 




V. 



(kj) 



(t) + H.c. 



(13) 



Here, we specify the initial state k and the field / of the per- 
turbation Vj.' and the unperturbed Hamiltonian Hq . 

The latter is defined as fff^ee i ^/ ~ ^k and contains cj^ as 
the frequency of level k G {0,1,2} and A/ as the detuning of 
field / e {1, 2}. We use a projector P^ — \k) {k\ on the levels 

k to identify coherent drive terms V"^' — V'--^^Pk starting 
from an initial state k. The superscript / e {1, 2} is used to 

split V{t) into VJ. for those terms which depend on ili and 

Vj. ' for the ones with il,2', a sign (±) denotes whether the 
initial process is an excitation (+), i.e. a term containing a 
factor e^'"^*, or a de-excitation (— ), with a factor e+*"/*. 

Using this formalism we find a considerable number of 
terms, time-independent and -dependent ones, some closer to 
resonance and others stronger detuned. Neglecting the time- 
varying terms rotating at twice a detuning A1/2 we obtain the 



effective two-photon Hamiltonian 

i = l,2/=l,2 *■ ^ ^' 



n} 



2{5i -62- A/) 



(|1>.(1|-|2),(2|) 
(|2),(0| + |0),(2|) 



4V2(5i -62- A 

Setting the detunings of the fields to A1/2 — Tlt^i + e) we 
have that Ai + A2 = and keep a certain virtual character of 
the single fields by a detuning of ±e, as can also be seen from 
Fig. [lib). In this configuration, there exists an effective two- 
photon drive where the first field (with ili) drives the lower 
transition |0) ^ |1) and the second field (with 5I2) drives the 
upper transition. Expressing the resulting effective Hamilto- 
nian in terms of the anharmonicity (using 61 ~ ^ — A) we 
obtain 



ii^eff ~ E 



^\ 



^l 



f^2 V 4e 4(2A + ^2 + e) 



(|0),(0|-|1>,(1|) 



VlI 



2(.2 + .)'2(2A + .)'(|l>^-<^|-|^>^-<2|) 



+ ^(-l)M|2),(0| + |0),(2|) 



with an effective two-photon Rabi frequency 



(15) 
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^2+6 2A 

2A82^{A-<i) 



e 2A + (^2 + e 



2^2 e{52 + (){2A + e){2A + 82 + ^)' 



(16) 



From here we see that for the case of zero anharmonicity 
A = 0, i.e. for harmonic transmons, no effective two photon 
drive is possible. For A ^ 0, however, there exists a possi- 
bility of driving the transition |0) O |2). Note that the re- 
maining diagonal terms in Eq. ( [T5| ) represent shifts which can 
be compensated by suitable (minor) detunings of the fields. 
Their effect on Eq. (fT6]l can be considered very small. In this 
case we can express Hd,cS as a single coherent coupling of 
the transition |0) O- |2), 



HacS — 



n 



eff 



J = l,2 



(-1FI2), 



H.C. 



(17) 



We have thus obtained the coupling constant ficff of the effec- 
tive two-photon coupling we introduced in Sec. Ill] With this 
result we can turn to the derivation of the effective Lindblad 
operators for the engineered decay mechanisms used for the 
preparation of the singlet state. 



B. Engineered decay processes and their effective Llndblad 
operators 

To model the effective, dissipative evolution we use the 
same effective formalism as in the previous section to derive 
the effective Lindblad operators ll44l 



^eff 



L, 



EE 



H, 



NH 



V^'^^^'it), 



(18) 



/ 



with the perturbative coherent excitation V{t) from an initial 
state fc by a field /, and a non-Hermitian Hamiltonian 



H, 



(kj) 
NH 



H, 



(kJl 



2^ 



Ak.f) 



J^r,i^n 



(19) 



with the perturbed Hamiltonian H^ ' defined previously. 
We focus on the effective resonator decay process activated 
by the two-photon drive Hcs and followed by decay of a res- 
onator excitation L^- With Hq = ^froc + -^cav ^(0 — ^cff 
(n <C 62), and Lm = ^k we arrive at an effective Lindblad 
operator 



TK 



/7h|S)(00| 



A^ 



>j)(S|, 



(20) 



with effective decay rates of k+ and k~ . This operator rep- 
resents the dissipative mechanism we engineer to rapidly pre- 
pare the singlet state |S) from |00). In addition, it includes 
the loss processes at rates of k~ from |S) into other states 
\(j)j) e {|11), |To), |T, 1), |00,2)}. Note that here we have 
ignored some less important terms as their effect on the popu- 
lation of the singlet is small. 

We calculate k^ of Eq. (20 1, using the driving from |00) 
to |So) as given by Eq. (17 1, with a matrix element of 



^S. The dynamics of the excited state |So) is described 



by the non-Hermitian Hamiltonian in Eq. ( 19 1 which cou- 
ples I So) to |S)|1) through the resonator interaction H^^^, 
forming a coupled subspace. For the non-Hermitian Hamil- 
tonian i?NH ' ^^ ^^^ subspace which contains |So) and 
is reached by excitation from |S) with the two-photon drive 

Hcs, we define i^NH.So = -^nh" with 

i^NH.So = ^2|So)(So| + (5i + Je)|S)|l)(l|(S| + 

+ V2g(|S)|l)(So|+H.c.). (21) 

In order to keep the notation compact, we have written the 
Hamiltonian in terms of the complex detunings Sj = Sj — ^ 

and Sc — Sc ~ ^ combining the energy with the imaginary 
line width of the levels. For the inverted operator we find 



H, 



-1 

NH,So 



= s. 



2,cff 



|So)(So| + <5i-:eff|S)|l)(l|(S|- 



-g-ffi(|S)|l)(So|+H.c.). 



(22) 



Here, we have introduced effective detunings of 62,cB ~ 



coupling constant of gcff — V^g ^ /i^ ° ■ Since the rate 

for resonator decay from | S) 1 1) into | S) is given by k, we gen- 
erally find an effective decay of k+ = ,- 'jl from |00) to 
|S), concluding that the effective coupling rate gcs governs 
the strength of the engineered decay process. 

The decay rate k+ is maximized by a parameter choice of 
62 ~ V^g and Sc = S2 — Si, which corresponds to the two- 
photon drive from |00) being in resonance with |So) and the 
resonator being resonant with the upper transitio n. W e then 



obtain ^eff 



and thus kj 



20„ 



In Sec. 







we will 



make use of this result to derive the error and the speed of the 
protocol. 

We now turn to the effective loss processes k~ as they ap- 



pear in Eq. (20 1. These rates can be calculated using the same 
procedure we applied for the derivation of k^ above. As |S) 
is driven to |Ti) by the two-photon drive we need to consider 

the non-Hermitian Hamiltonian i^NH.Ti = ^nh "" which 
describes the subspace consisting of |Ti) and the states cou- 
pled to it by -fff-av Foi" low anharmonicities A < S2, i^NH.Ti 
needs to reflect the full complexity of the coupled subspace 
containing |Ti), |11>|1), |To)|l), |T)|2) and |00)|3). For an- 
harmonicities of j4 > $2, however, the subspace of |Ti) and 
1 11) 1 1) begins to decouple from the other states so that the dy- 
namics of the excited states can be approximated using only 
|Ti) and |11)|1). The Lindblad operator of Eq. ^ for the 
effective resonator decay then reduces to 



^off 



IS) (001 



11) 



(23) 



containing a single loss rate k_ = k^j^, from |S) into |11). 

To derive k_, we then approximate i^NH.Ti by the non- 
Hermitian Hamiltonian of the excited subspace consisting of 

|Ti)and|ll)|l), 

HnH.T, « ^2|Ti)(Ti| + (^1 + Je)|ll)|l)(l|(ll| + 

+ 25(|ll)|l)(Ti|+H.c.), (24) 

using the complex detunings defined above. The inverted op- 
erator is then given by 



H, 



NH,Ti 



4:eff|Ti)(Tii + ^rc!cff|ii>|i>(i|(T|+ 

+ 52^0(1 (|ll)|l)(Ti|+H.c.). (25) 



Here, we have found effective detunings 62.cii.T1 = ^2 — j- 



and (5ic,cff,Ti = (^1 + Sc) 



Si+S, 



and an effective coupling 



constant of geS.Ti —'^9 2 ' which are different from 

the ones in the previous case of |So). With a decay rate k from 
|11)|1) into 1 11), we obtain an effective decay rate of k_ « 

"" for the losses from |S). For the above choice of 82 

9 
which results in k_ w '^"'-o" . From here we conclude that 



ISott.TiP 

and (5c, the effective coupling constant becomes .gcff.Ti 

4ff2 

for K^ <C g^ the effective loss rate k_ from the singlet is 
engineered to be much smaller than its preparation rate k+ « 



<52 - ^ and (5ic,off = (^1 + 4) - j^, and an effective ?^. These results confirm the explanations in Sec. 
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FIG. 3: (Color). Dressed state energy vs. anharmonicity. An ef- 
fective two-photon drive f^cft from |00) (solid black line) to jS_) 
(green dashed) is implemented as two consecutive single-photon ex- 
citations by two microwave fields, Q± and fl2- The individual drives 
are mediated by a dressed state of |T) (short-dashed red) and made 
virtual through a detuning of e (not shown). The two virtual exci- 
tations combine to an effective drive Qcft resonant with the transi- 
tion |00) -f-^ jS-); jS-) then decays into |S) (indicated). The same 
field couples to the transition from |S-) to the dressed states of |Ti) 
(dashed-dotted). By an appropriate choice of the oscillator detuning 
(5c (here plotted for Sc — S2 — <5i with a; — 20g), this coupling to 
|Ti) is made off-resonant (left set of arrows). In case that |Ti) is hit 
by the drive (right set of arrows), 5c is chosen differently to make the 
coupling off-resonant. 



Note that, on the one hand, the above treatment of the cou- 
pled excited subspace where we restrict the excited state sub- 
space to |Ti) and |11)|1) is quite simplistic, given that it re- 
duces the number of resonances from five to only two. In 
particular, one needs to ensure that one does not hit an acci- 
dental resonance with one of the dressed states of the system. 
On the other hand, the parameter space consisting of Sc, 82 
and e is sufficiently big to avoid an excitation of the remain- 
ing undesired resonances as there are sufficiently many suit- 
able points in different regions of parameter space for which 
all of these resonances are off-resonant with the two-photon 
drive. In Fig. l3] we draw the dressed states of the coupled 
resonator-transmon system. Here, the single-photon fields are 
tuned to resonantly excite the transition |00) O |S_) by a 
two-photon transition, mediated by the triplet state |T). The 
same two-photon drive also couples | S) to a number of dressed 
states with contributions from iTi). These transitions, how- 



ever, generally have different frequencies than the desired 
one from |00) to |S_) so that excitation of |S) by the drive 
ricff is off-resonant and suppressed by its detuning from the 
dressed states. This can be seen from Fig. [3] where we draw 
the dressed states together with the two-photon drive for the 
choice of 5c = 62 — 5i. In the figure, we show an example 
near A ^ '^ where the driving is off-resonant with the excited 
states which contain contributions from |Ti). We also draw 
an example at A f^ 2g where this is not the case and where a 
resonance is hit accidentally. Here, it is necessary to choose a 
different detuning Sc- Below, we will verify by numerical sim- 
ulation for a broad parameter range that it is always possible 
to avoid such resonances. 

In addition to losses caused by the two-photon drive, also 
the individual fields ili and 0,2 couple |S) to other states. The 
coupling of the even-phase single-photon drive ili from |00) 
to |So) does not cause any significant loss from |S), since pop- 
ulation in |So) is recycled via |S)|1) back into |S). The odd- 
phase single-photon drive ^2, on the other hand, couples |S) 
to |00) and to a superposition state ^TjdH) - |To)). Both 
these states are dark states of the resonator coupling. Thus, no 
exchange excitation to a resonator-excited state can shift them 
into resonance with the off-resonant drive $12 from |S) and no 
effective resonator decay process from |S) is established in- 
volving them. Accumulation in these states does not occur, 
either, given that 4=(|11) — |To)) decays through qubit de- 
coherence and 1 00) decays into |S) as discussed earlier. As a 
consequence, neither of the two drives causes significant loss 
from the singlet. 

Another source of errors emerges for small anharmonici- 
ties A < S2 from the coherent coupling of |S) to other states 
like 1 00) by the single-photon drives Oi and il2- However, 
for A > 62, these couplings are sufficiently detuned to be ig- 
nored. Also, beside effective resonator decay processes, qubit 
decoherence occurs according to Eqs. ([8])-(|9]l. Provided that 
the decay rate 7 is much weaker than all other physical cou- 
plings present in the system, i.e. j -^ K,g, effective processes 
combining qubit decoherence with coherent excitation can be 
safely neglected. 

We conclude that the sources of error originating from ef- 
fective resonator decay which can cause losses from the sin- 
glet state are suppressed for the right parameter choice. These 
processes are, together with the engineered dissipative state 
preparation process, contained in the effective resonator de- 



cay operator in Eq. ( 23 1 



IV. PARAMETER AND PERFORMANCE ANALYSIS, 
IMPERFECTIONS AND REALIZATION ASPECTS 

In the previous section we have identified the effective co- 
herent and dissipative processes which are relevant for our 
dissipative state preparation scheme and investigated the cor- 
responding Lindblad operators and rates. In this section, we 
will use these results to derive approximate expressions for the 
error and speed of the presented protocol as the main bench- 
marks for our scheme. Later, we will assess the temporal evo- 
lution of the system numerically. 



A. Error and speed of the protocol 

In the previous section we have derived the effective res- 
onator decay operator L^jj, given in Eq. (23 1, which describes 



both the preparation of and the losses from the singlet state 
|S). The derivation of Eq. (23 i was carried out in the limit of 



weak driving. As we will find numerically below, the dissipa- 
tive preparation of the singlet at a rate of k+ « 20^^^/^ works 
well for a driving strength up to ftcs ~ K./4, which yields a 
preparation rate k+ « k/8 for the singlet state |S) and a loss 
rate k_ « k"^ /6Ag^ from it. In addition, |S) decays at a rate 
of 7, as described by the operators in Eq. (|8]l-(j9]l. 

Based on these rates we can approximate the temporal dy- 
namics with rate equations of the populations Pi = {ipi\p\ijji) . 
We assume that the reshuffling mechanism rapidly transfers 
all population from the triplet states to the state |00), which is 
correct for small anharmonicity A, the evolution of the pop- 
ulation of the singlet can then be described by a single rate 
equation for the population of the singlet Ps, 



Ps 



K+Pm 



(k_+7)Ps, 



(26) 



formulated in terms of the decay rates specified above. By 
simply comparing the gain and loss of the singlet in the steady 
state, i.e. Pg = 0, we can estimate the steady-state fidelity 
Fs = lim Pg of the singlet and, consequently, the error of 

the protocol (1 — Pg)- Assuming a near unit fidelity we obtain 



(1-Ps) 



7 + K- 



87 

K 



(27) 



From this expression we can readily see that the error of the 
protocol has a promising scaling with the physical parameters. 
Specifically, the error depends on the ratios of coupling and 
noise, g/n and K/7 so that it will be small for strong coupling, 
g^ ^ K^, and modest qubit decoherence, 7 -^^ k. Under the 
assumption that we can vary the resonator decay rate k we can 
minimize the error in Eq. ( [27] i by choosing k. Considering 
g^(l — Pg) = 0, we derive the optimal resonator decay rate 

2 ^/4^yg^. Inserting this yields the optimized error of 



''opt 



the protocol. 



(1 - Ps) 



opt 



2/3 



(28) 



From here we conclude that for 7 <^ g the error of the proto- 
col can be limited to very small values at the few percent level 
or even lower for realistic parameters. We will later confirm 
this finding numerically. 

In addition, the convergence time, i.e. the decay time of the 
undesired states, can be approximated using Eq. (J26]l, assum- 
ing rapid reshuffling of the undesired states to 1 00) . Given that 
here the preparation of the singlet at a rate k+ is the dominant 
process, the convergence time r is given by 
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FIG. 4: (Color) Evolution of the system towards an entangled steady 
state. Initially prepared in an equal mixture of the lower states (|00) - 
green, dotted line, 1 11) - red, dashed-dotted line, |T) - blue, dashed 
line, |S) - purple, solid line) the system evolves towards its steady 
state which is close to the maximally entangled singlet state of the 
two transmons. Part a) and b) show the result for an anharmonicity of 
A = g and A = 4.75(/ respectively. The inset in a) shows the region 
in the Aa — Ag plane where the singlet fidelity is high, Fs > 90%. 
For A = g the number on each contour line indicates the preparation 
time in units of 1/g. The inset in b) shows the singlet state fidelity 
at i = 1000/(7 2S a function of anharmonicity. Parameter values are 
^1/2 = 5/3, fi: ~ Sg/lO and 7 — ,g/5400. The values of lj, A1/2 
and Sc are obtained through numerical optimization. 



resonator decay rate k. 

Note that the above expressions for the error and the con- 
vergence time are approximative and derived using results ob- 
tained for the assumption of weak driving in Sec. IIIB In 



our numerical simulations below we will optimize a number 
of parameters including the driving strength to achieve highly 
entangled states within a preparation time as short as possible. 
In doing so, we arrive at particular choices of the available 
parameters which allow us to achieve high fidelities in short 
time. As these optimal parameters are in a regime where the 
effective Lindblad operators do no longer accurately describe 
the dynamics |44|, the findings of Eqs. (27i and (29i deviate 
from the simulation results below. 



(29) 



B. Numerical results 



where we have used iloff ~ '«/4. We thus conclude that the 
timescale of dissipative state preparation is determined by the 



To verify the findings above as well as to investigate the 
limitations of the approximation we now depart from the an- 
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alytical treatment in the previous sections and assess the per- 
formance of the scheme numerically |45|. To this end we inte- 
grate the master equation in Eq. dTJi including the four lowest 
levels of each transmon and up to three photons in the res- 
onator Due to the Stark shifts induced by the driving, we have 
numerically optimized the sum- and difference frequencies w 
and A 1/2 of the drives, as well as the resonator frequency Sc- 
In Fig. Blwe plot the populations 



P,(<)=Tr((|l^,;)(V',|®lcav)p(i)) 



(30) 



between the time evolved density matrix p{t) and the four 
lower states \ipi) = |00), |11), |S), |T) inti'oduced in Sec. [S] 
The results of our simulation are shown in Fig. |4]a)-b), where 
we plot the populations, starting with an initial, equal mixture 
of all four lower states. In Fig. l4]a), we consider a rather low 
anharmonicity A = g, which is also what is typically used 
in experiments 191 (TT] [12 1 . Here, the population of the states 
|11) and |T) show a fast drop due to the reshuffling into |00). 
At the same time, albeit on a slightly longer timescale, the 
dissipative preparation of the singlet is performed, reaching a 
fidelity of 90% within a time of about r « "200/ g. For a trans- 
mon experiment with g/{2TT) — 300 MHz this would allow 
preparation times of about r sa 80 ns. For the results in Fig. 
[4] we have chosen 7/(27r) sa 60 kHz w (7/(27r5400) corre- 
sponding to a relaxation time of Ti w 3 /xs L2J for the above 
parameter choice. This is much shorter than current state-of- 
the-art transmon qubits which have decoherence times of up 
to T2 ^ 95/is and Ti ^ TQjjls lH |9l and thus very conserva- 
tive. We chose this somewhat low number for the transmon 
relaxation time to show the robustness of our scheme against 
such imperfections and to demonstrate the possibility of im- 
plementing it not only in state-of-the-art 3D, but also in 2D 
transmon systems. We find that, with state of the art numbers 
for transmon decoherence times as quoted above, it is possible 
to reach a steady state fidelity of ~ 97% for A — g. 



C. Anharmonicity of the transmon 

As discussed in the previous sections, the coupling of the 
resonator to the |0) -O- |1) transition for each transmon con- 
tributes to the scheme by reshuffling the unwanted populations 
to 1 00). This coupling, however, gets increasingly detuned for 
higher anharmonicities A. In Fig. |4]b) we show the effect of 
an increasing A on the preparation scheme. Here, for a rather 
high anharmonicity of A = A.lbg, the reshuffling of the states 
|11) and |T) to |00) is slowed down as compared to the result 
for A = g in Fig. |4]a). This can be seen from the drop in the 
population of |T) and |11) which is much less pronounced in 
b) than in a). In addition, we observe an increase in the steady 
state populations of these states. It is therefore advantageous 
to work with a rather low anharmonicity, where the coupling 
to the lower transition is still effective. Such anharmonicities 
are typical for state-of-the-art experiments ll9l [m[T2l . 

In the following, we will assess the possibility to operate 
our scheme for a broader range of anharmonicities, despite 
the breakdown of the reshuffling. To this end we allow for a 
rather long preparation time t = IOOO/5. In the inset in Fig. 



|4]b) we show results achieved using a numerical optimization 
routine to optimize the fidelity by fine-tuning the frequencies 
of the microwave fields and the resonator. These degrees of 
freedom in the parameter choice are used by the optimization 
routine to avoid undesired resonances by a slight departure 
from the resonance conditions of the previous sections. The 
range of our protocol is then limited by the breakdown of the 
reshuffling to A < 4g, as well as to A > g. For lower A the 
effective two-photon drive becomes ineffective and couplings 
to higher levels of the transmons add shifts to the resonances 
required for the state preparation mechanism. To reach a high 
fidelity Fs > 90% of the steady state one should therefore 
work with anharmonicities between A ^ g and A « Ag. 

Finally, we briefly comment on the possibility for dissipa- 
tive state preparation with even more anharmonic systems: In 
this case we choose to have the resonator in (or close to) res- 
onance with the upper transition. Consequently, the lower 
transition is largely detuned and its effect negligible. We 
thereby achieve a situation which is very similar to optical 
resonator QED with atomic A schemes - a system where var- 
ious schemes for dissipative preparation of entanglement are 
available Il23ll25l . These schemes can then be mapped to the 
highly anharmonic circuit QED setup. In those schemes the 
role of the far-detuned resonator coupling on the lower transi- 
tion is accomplished by an additional microwave field which 
takes over the reshuffling of the triplet states. In this way, 
preparation of a steady state close to the maximally entangled 
singlet state can be achieved for any anharmonicity. For low 
anharmonicitiy, however, the coupling of the resonator to the 
lower transition allows us to avoid this field and thus to sim- 
plify the experimental implementation. 



D. Experimental imperfections 

From the previous discussion it is clear that our scheme re- 
lies on the fact that the two transition frequencies of the trans- 
mons are identical. Moreover, we have so far only considered 
the case when the coupling, g, is identical for both transmons. 
In this section, we depart from these assumptions and consider 
the effect of experimental imperfections. The transmons are 
characterized by their spectrum which is is set by the effective 
Josephson energy, Ej and the charging energy Ec = 2 A BT| . 
Since Ej can be tuned in situ we assume that the lower transi- 
tion frequencies are tuned into resonance. We therefore focus 
our analysis on the charging energy (anharmonicity) and the 
couplings by considering A2 — A^^i and 32 — ^g9i where 
the subscript denotes transmon number. In the inset of Fig. 
pk, we plot the region in the A^i — Ag plane where F > 90% 
for A = g. The different contours correspond to the indicated 
preparation time and we see that there is roughly a 20% error 
tolerance built into the system with respect to these parame- 
ters. The reproducibility of Ec and g is set by the precision of 
the e-beam lithography process and these tolerances are well 
within the limits of current technology. 
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V. CONCLUSION AND OUTLOOK 

In this work we have presented a scheme for the prepara- 
tion of an entangled steady state of two transmons by means 
of dissipation. We have engineered effective decay mecha- 
nisms for the dissipative preparation of the desired maximally 
entangled singlet state and verified them analytically and nu- 
merically. We have demonstrated that high fidelity with the 
singlet state can be reached within favorable time for realistic 
experimental parameters, even in the presence of rather strong 
decoherence of the qubits. In addition, our scheme has proven 
to be robust against experimental imperfections such as non- 
degeneracy of the transmon levels and couplings. We hope 
that our proposal will be an important stepping stone for the 
generation of high-fidelity and stationary entanglement gener- 
ation in circuit QED setups and for more sophisticated dissi- 
pative protocols in superconducting systems. 
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Note added. Under the final stage of manuscript prepara- 
tion, our attention was drawn to the submission of a study 
|!46l with a similar objective. A full comparison to this work 
is beyond the scope of our paper. We note, however, that the 
coherence times assumed there are much longer than in our 
work. Furthermore, due to the highly dispersive nature of the 
interaction assumed in that study, the convergence time pre- 
sented there is significantly longer. 
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